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Abstract 

We study Radiative Electroweak Symmetry Breaking in the Minimal Supersymmetric Standard 
Model (MSSM). We employ the 2-loop Renormalization Group equations for running masses 
and couplings taking into account sparticle threshold effects. The decoupling of each parti- 
cle below its threshold is realized by a step function in all one-loop Renormalization Group 
equations (RGE). This program requires the calculation of all wavefunction, vertex and mass 
renormalizations for all particles involved. Adapting our numerical routines to take care of the 
succesive decoupling of each particle below its threshold, we compute the mass spectrum of 
sparticles and Higgses consistent with the existing experimental constraints. The effect of the 
threshold corrections is in general of the same order of magnitude as the two-loop contributions 
with the exception of the heavy Higgses. 
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The purpose of the present talk is to report briefly on a treatment of low energy threshold 
effects in the Renormalization Group equations of the parameters of the MSSM in the framework 
of Radiative Electroweak breaking. Since we have employed the DR scheme in writting down 
the one- loop Renormalization Group equations, which is by definition mass-indepedent, we 
could "run" them from Mx down to Mz without taking notice of the numerous sparticle 
thresholds existing in the neighborhood of the supersymmetry breaking scale near and above 
Mz- This approach of working in the "full" theory consisting of particles with masses varying 
over 1-2 orders of magnitude has to overcome the technical problems of the determination of 
the pole masses. Our approach, also shared by other analyses, is to introduce a succession of 
effective theories defined as the theories resulting after we functionally integrate out all heavy 
degrees of freedom at each particle threshold. Above and below each physical threshold we 
write down the Renormalization Group equations in the DR scheme only with the degrees 
of freedom that are light in each case. This is realized by the use of a theta function at 
each physical threshold. The integration of the Renormalization Group equations in the "step 
approximation" keeps the logarithms ln(— ) and neglects constant terms. The physical masses 
are determined by the condition m(m p h ys ) = ni v hys which coincides with the pole condition if 
we keep leading logarithms and neglect constant terms. The great advantage of this approach 
is that the last step of determining the physical mass presents no extra technical problem and 
it is trivially incorporated in the integration of the Renormalization Group equations. 

A dramatic simplification of the structure of the supersymmetry breaking interactions is 
provided either by Grand Unification assumptions or by Superstrings. The simplest possible 
choice at tree level is to take all sparticle and Higgs masses equal to a common mass parameter 
m , all gaugino masses equal to some parameter m 1( / 2 and all cubic couplings flavour blind and 
equal to A . This situation is common in the effective Supergravity theories resulting from 
Superstrings but there exist more complicated alternatives. For example Superstrings with 
massless string modes of different modular weights lead to different sparticle masses at tree 
level@. The equality of gaugino masses can also be circumvented in an effective supergravity 
theory with a suitable non-minimal gauge kinetic termB. Note however that such non- minimal 
alternatives like flavour dependent sparticle masses are constrained by limits on FCNC pro- 
cesses. In what follows we shall consider this simplest case of four parameters m , m^, A D 
and B . The scalar potential of the model is 

V = mllH^ 2 + m\\H 2 \ 2 + # 2 + c.c) 

+ l -g' 2 {\H x \ 2 - \H 2 \ 2 ) 2 + + |tf 2 | 4 + A\HlH 2 \ 2 - m^Hil 2 ) + (1) 

o o 
written in terms of 

m h = m H lt2 + ( 2 ) 

Only an unbroken minimum appears for m\ = m\. Replacing the appearing parameters with 
their running values m 2 (Q), m^Q),... as defined by the Renormalization Group and adding 
the one-loop radiative corrections obtained in the DR scheme, 

AVi = ^- 2 Str{M\ln(M 2 /Q 2 ) - 3/2)} (3) 

we end up with an Effective Potential that upon minimization supports a vacuum with spon- 
taneously broken electroweak symmetry^™. A reasonable approximation to (3) would be to 
allow only for the dominant top-stop loops. Note that although the Renormalization Group 
improved tree level potential depends on the scale Q this is not the case for the full 1-loop 
Effective Potential which is Q-independent up to, irrelevant for minimization, Q-dependent but 
field-independent terms. 



We shall assume that at a very high energy scale Mx the soft supersymmetry breaking 
is represented by four parameters m , mi/2, A Q and B of which we shall consider as input 
parameters only the first three and treat B(Mz) as determined by minimization conditions of 
the one loop effective potential. Actually we can treat (3(Mz) as input parameter and both 
B(Mz), n(Mz) are determined by solving the minimization conditions with the sign of \x left 
undetermined. The top-quark mass^t, or equivalently the top-quark Yukawa coupling, although 
localized in a small range of values should also be considered as an input parameter since the 
sparticle spectrum and the occurrance of symmetry breaking itself is sensitive to its value. Thus, 
the input parameters are m , m 1 / 2 , A Q , f3(M z ) and m t (M z ) as well as the sign of \i. 

In our notation, for a physical mass M, 

9 M = 9(Q 2 - M 2 ) (4) 

Also t stands for t = InQ 2 and (3\ = % for each parameter A. Note also that we assume 
diagonal couplings in family space. 

As an example consider the one loop /3-function of the trilinear couplingll A T , 

^ = j^{-3g 2 2 M 2 9 Wfli -^g 1 2 M 1 (2 + 9 f{l )9 s 

+ 3r b %^Q + 4F r 2 A T + A T [Z rl ^ 2 + Z T2 ^ 2 + Z TT r T 2 ]} (5) 

Where Z rl , Z r2 and Z TT are displayed in Table I. 

TABLE I 

Z Tl = A [11 + 10^ - %9 & - A9^ + 89^ + 

26 Hl - 89 Hi e - 26 L - 6g L - 89 %i - ^b Hi L + ^h.l] 
Z t2 = l[-3 + 69 Hl - 69 L - 129 HiL + 69 w - 39 LW + I29 wilii ] 
Z TT = ±[-16 + +69 &i - 9^ - 39 Hl + 49 Hi z + 4B 6L - 29 kll + 89 Hii ] 



Table I: Threshold coefficients appearing in the renormalization group equation 
of the trilinear scalar coupling A T . 

Note that the threshold corrections introduced in our approximation by the theta-functions 
at 1-loop are expected to be comparable to the standard 2-loop RG corrections. In our nu- 
merical analysis that we follow we shall employ the 2-loop RG equations which have not been 
presented here due to their complicated form but can be found elsewhere 0. The problem at 
hand consists in finding the physical masses of the presently unobserved particles, i.e. squarks, 
sleptons, Higgses, Higgsinos and gauginos, as well as their physical couplings to other observed 
particles. This will be achieved by integrating the Renormalization Group equations from a 
superheavy scale Mx, taken to be in the neighbourhood of 10 16 GeV, down to a scale Q in 
the stepwise manner stated. If the equation at hand is the Renormalization Group equation 



for a particular running mass m(Q), then Q Q is the corresponding physical mass determined 
by the condition m(Q ) = Q D . If the equation at hand is the Renormalization Group equa- 
tion for a coupling the integration will be continued down to Q D = Mz- Acceptable solutions 
should satisfy the minimization conditions at Mz, i.e. describe a low energy theory with broken 
electroweak symmetry at the right value of Mz — 91.187CreV. 

The boundary condition at high energy will be chosen as simple as possible, postponing 
for elsewhere the study of more complicated alternatives. Thus at the (unification) point Mx, 
taken to be 10 16 GeV, we shall take 

m Q {M x ) = m 5c (M x )=m Dc (M x )=m L (M x )=mEc(M x ) 

= m Hl (M x ) = m H2 (M x ) = m (6) 

and 

M 1 (M X ) = M 2 (M X ) = M 3 (M X ) = mi/a (7) 
In addition we take equal cubic couplings at M x , i.e. 

A e (M x ) = A d (M x ) = A U (M X ) = A (8) 

Our set of constraints includes the low energy experimental gauge coupling values which 
we have taken to be M z = 91.1876?eu , a(M z )~ 1 jYs = 127 - 9 ± O- 1 and ( sin2 °w)ws = 0.2316 - 
.881(T 7 (M t 2 - 160 2 )Gev- 2 . The knowing, average experimental value of a 3 is 0.117 ± 0.010. 
These MS values for the couplings are related to the relation, DB^pnes through the relations 
9jTs = 9Tm(^ ~ Cg 2 /96n 2 ), where C = 0, 2, 3 respectively for the three factor gauge groups. For 
the b-quark and r-lepton masses we have taken m& = 5.0 Gev and m T = 1.8 Gev.The recent 
evidence!!' for the top quark mass has motivated values in the neighborhood of 176 ±8 Gev.The 
physical top quark mass M t is related to the running top-quark mass through the approximate 
relation 

= (i + g + ,„.) < 9 > 

As stated previously the B, \x are not inputs in the approach we are following but are 
determined through the equations minimizing the scalar potential. For their determination at 
the scale Mz we take into account the one loop corrected potential considering the dominant 
top and stop contributions. This procedure modifies the tree level values B(M Z ), fJ-(M z ) ■ It is 
well known that the value of \x affects the predictions for the physical masses especially those 
of the neutralinos and charginos. In approaches in which the effect of the thresholds is ignored 
in the RGE's the determination of B,li is greatly facilitated by the near decoupling of these 
parameters from the rest of the RGE's. However with the effects of the thresholds taken into 
account such a decoupling no longer holds since the thresholds themselves depend on B,fi, 
or equivalently on //, m|. Thus, as initial inputs for B(M Z ) and \i{Mz) we take those arising 
from the minimization equations assuming that theshold effects are absent. At this stage our 
analysis is identical to those of other authors. Subsequently we run our numerical routines 
switcing on the threshold contributions to the RGE's keeping fixed the inputs for A ,m , mi/2, 
tan/? and all couplings. This procedure corrects the initial inputs for B(Mz), \i{Mz) in each 
run until convergence is reached. This is unecessary of course in cases where the thresholds are 
neglected. The next step regarding the mixing parameters /a, m\ is to correct them taking into 
account the one loop effective potential in the way prescribed earlier. 



1 Note that at the 2-loop order the DR scheme needs to be modified so that no contribution to the scalar 
masses due to the"e-scalars" lEl shows up. 



We have displayed some of our results in tables II and III. We have taken n(M z ) positive 
and hence BiMz) negative. Their mirror values \i{M%) < 0, BiMz) > lead to qualitatively 
similar results. In the table II, for a characteristic set of values A Q = 400 GeV, m = 300 
GeV and = 200 GeV we have varied tan(3 between 2 and 25. Note the well known© 
approximate equality between the masses of one of the neutralinos and one of the charginos. 
The lightest Higgs turns out to be heavier than the Z - boson. Althought not displayed, for 
negative fi its mass drops below Mz for small values of the angle tanf3 ~ 2. 

Finally, Table III compares for a characterestic choice of parameter values, three dinstict 
cases. Case [a] indicates one loop predictions, case [b] two loop predictions with thresholds only 
in couplings, and case [c] the complete two loop with thresholds everywhere. Comparison of the 
first two cases [a] and [b] point out the fact that thresholds in couplings affect only by small 
ammount (1% — 2%) the spectra, except from the neutralino and chargino states, where the 
differences are of order 10% due to different evolutions of the soft gaugino masses M\^- 

Comparing cases [b] and [c], we observe quite large effects in states labeled as X34, Xi 
as well as in heavy Higgs states. This discrepancy is due mainly to the evolution of m|, whose 
values affect substantially the masses of the pseudoscalar and charged Higgses,and in particular 
on it's dependense on the gaugino masses ILil . 



TABLE II 

m t = 175 , A Q = 400 , m = 300 , m 1/2 = 200 , fi(M z ) > 
tanp 25 20 15 10 2 



M GUT 2.632 2.633 2.633 2.632 2.515 
(10 16 GeV) 

a GUT .04161 .04161 .04161 .04161 .04134 

127.9 127.9 127.9 127.9 127.9 

sin 2 9 w .2311 .2311 .2311 .2311 .2311 

a 3 .13128 .13129 .13129 .13126 .12909 



M t 177.0 177.0 177.0 177.0 176.8 



g 495.6 495.8 495.9 495.9 492.4 

Xi 77.1 77.0 76.8 76.4 74.7 

X° 2 139.1 138.9 138.5 137.6 136.0 

X° 3 350.5 352.3 354.7 359.1 487.8 

X° A -337.2 -338.6 -340.3 -343.2 -467.7 

X\ 352.3 354.0 356.0 359.9 484.5 

X| 13^8 138^3 138J 137X) 134.7 



t u t 2 527.3,315.4 531.9,315.2 535.9,314.4 539.4,312.3 543.9,285.3 

61,62 506.2,441.0 514.2,451.4 520.9,459.4 525.8,465.0 525.9,462.4 

fi,f 2 328.0,266.5 330.4,281.0 331.7,292.9 331.9,302.1 329.4,309.2 
z> r 306.4 311.5 315.6 318.6 323.4 



ui,2,«i >2 537.8,528.9 537.8,528.9 537.8,528.9 537.7,528.8 535.4,525.8 
^1,2,^1,2 543.3,529.6 543.3,529.6 543.3,529.6 543.2,529.5 538.6,525.8 

e 1>2 ,e c 1>2 330.1,311.7 330.1,311.7 330.0,311.7 330.0,311.6 328.8,310.3 
i/1,2 ' 321.0 321.0 320.9 321.0 323.5 



A 630.7 612.0 588.1 560.1 669.0 

h ,H 114.1,630.6 114.2,611.9 114.2,588.0 113.8,560.2 92.1,672.8 
H± 635.4 616.9 593.1 565.4 673.5 



Table II: MSSM predictions for m t = 175 GeV, A = 400 GeV, m = 300, m 1/2 = 200GeV 
and for values of tan (3 ranging from 2 to 25. Only the \x > case is displayed. 



TABLE III 
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Table III: MSSM mass spectrum for the inputs shown in the first row (// > 0). We compare 1 
- loop (case [a]), 2 - loop with thresholds in couplings (case [b]) and complete 2 - loop predictions 
(case [c]) with thresholds in both couplings and dimensionful parameters. 
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